Interpretation of morphogen gradients by a synthetic bistable circuit.

During development, cells gain positional information through the interpretation of dynamic morphogen gradients. A proposed mechanism for interpreting opposing morphogen gradients is mutual inhibition of downstream transcription factors, but isolating the role of this specific motif within a natural network remains a challenge. Here, we engineer a synthetic morphogen-induced mutual inhibition circuit in E. coli populations and show that mutual inhibition alone is sufficient to produce stable domains of gene expression in response to dynamic morphogen gradients, provided the spatial average of the morphogens falls within the region of bistability at the single cell level. When we add sender devices, the resulting patterning circuit produces theoretically predicted self-organised gene expression domains in response to a single gradient. We develop computational models of our synthetic circuits parameterised to timecourse fluorescence data, providing both a theoretical and experimental framework for engineering morphogen-induced spatial patterning in cell populations.

Turing Patterning in Stratified Domains.

Reaction-diffusion processes across layered media arise in several scientific domains such as pattern-forming E. coli on agar substrates, epidermal-mesenchymal coupling in development, and symmetry-breaking in cell polarization. We develop a modeling framework for bilayer reaction-diffusion systems and relate it to a range of existing models. We derive conditions for diffusion-driven instability of a spatially homogeneous equilibrium analogous to the classical conditions for a Turing instability in the simplest nontrivial setting where one domain has a standard reaction-diffusion system, and the other permits only diffusion. Due to the transverse coupling between these two regions, standard techniques for computing eigenfunctions of the Laplacian cannot be applied, and so we propose an alternative method to compute the dispersion relation directly. We compare instability conditions with full numerical simulations to demonstrate impacts of the geometry and coupling parameters on patterning, and explore various experimentally relevant asymptotic regimes. In the regime where the first domain is suitably thin, we recover a simple modulation of the standard Turing conditions, and find that often the broad impact of the diffusion-only domain is to reduce the ability of the system to form patterns. We also demonstrate complex impacts of this coupling on pattern formation. For instance, we exhibit non-monotonicity of pattern-forming instabilities with respect to geometric and coupling parameters, and highlight an instability from a nontrivial interaction between kinetics in one domain and diffusion in the other. These results are valuable for informing design choices in applications such as synthetic engineering of Turing patterns, but also for understanding the role of stratified media in modulating pattern-forming processes in developmental biology and beyond.

Programmed hierarchical patterning of bacterial populations.

Modern genetic tools allow the dissection and emulation of fundamental mechanisms shaping morphogenesis in multicellular organisms. Several synthetic genetic circuits for control of multicellular patterning have been reported to date. However, hierarchical induction of gene expression domains has received little attention from synthetic biologists, despite its importance in biological self-organization. Here we report a synthetic genetic system implementing population-based AND-logic for programmed autonomous induction of bacterial gene expression domains. We develop a ratiometric assay for bacteriophage T7 RNA polymerase activity and use it to systematically characterize different intact and split enzyme variants. We then utilize the best-performing variant to build a three-color patterning system responsive to two different homoserine lactones. We validate the AND gate-like behavior of this system both in cell suspension and in surface culture. Finally, we use the synthetic circuit in a membrane-based spatial assay to demonstrate programmed hierarchical patterning of gene expression across bacterial populations.

Orthogonal intercellular signaling for programmed spatial behavior.

Bidirectional intercellular signaling is an essential feature of multicellular organisms, and the engineering of complex biological systems will require multiple pathways for intercellular signaling with minimal crosstalk. Natural quorum-sensing systems provide components for cell communication, but their use is often constrained by signal crosstalk. We have established new orthogonal systems for cell-cell communication using acyl homoserine lactone signaling systems. Quantitative measurements in contexts of differing receiver protein expression allowed us to separate different types of crosstalk between 3-oxo-C6- and 3-oxo-C12-homoserine lactones, cognate receiver proteins, and DNA promoters. Mutating promoter sequences minimized interactions with heterologous receiver proteins. We used experimental data to parameterize a computational model for signal crosstalk and to estimate the effect of receiver protein levels on signal crosstalk. We used this model to predict optimal expression levels for receiver proteins, to create an effective two-channel cell communication device. Establishment of a novel spatial assay allowed measurement of interactions between geometrically constrained cell populations via these diffusible signals. We built relay devices capable of long-range signal propagation mediated by cycles of signal induction, communication and response by discrete cell populations. This work demonstrates the ability to systematically reduce crosstalk within intercellular signaling systems and to use these systems to engineer complex spatiotemporal patterning in cell populations.

The identification of QTL controlling ergot sclerotia size in hexaploid wheat implicates a role for the Rht dwarfing alleles.

KEY MESSAGE: Four QTL conferring resistance to ergot were identified in the UK winter wheat varieties 'Robigus' and 'Solstice'. Two QTL co-located with semi-dwarfing alleles at the Rht loci Rht - 1B and Rht - 1D implicating a role of these DELLA proteins in infection success of Claviceps purpurea. The fungal pathogen Claviceps purpurea infects ovaries of a broad range of temperate grasses and cereals, including hexaploid wheat, causing a disease commonly known as ergot. Sclerotia produced in place of seed carry a cocktail of harmful alkaloid compounds that result in a range of symptoms in humans and animals, causing ergotism. Following a field assessment of C. purpurea infection in winter wheat, two varieties 'Robigus' and 'Solstice' were selected which consistently produced the largest differential effect on ergot sclerotia weights. They were crossed to produce a doubled haploid mapping population, and a marker map, consisting of 714 genetic loci and a total length of 2895 cM was produced. Four ergot reducing QTL were identified using both sclerotia weight and size as phenotypic parameters; QCp.niab.2A and QCp.niab.4B being detected in the wheat variety 'Robigus', and QCp.niab.6A and QCp.niab.4D in the variety 'Solstice'. The ergot resistance QTL QCp.niab.4B and QCp.niab.4D peaks mapped to the same markers as the known reduced height (Rht) loci on chromosomes 4B and 4D, Rht-B1 and Rht-D1, respectively. In both cases, the reduction in sclerotia weight and size was associated with the semi-dwarfing alleles, Rht-B1b from 'Robigus' and Rht-D1b from 'Solstice'. Two-dimensional, two-QTL scans identified significant additive interactions between QTL QCp.niab.4B and QCp.niab.4D, and between QCp.niab.2A and QCp.niab.4B when looking at sclerotia size, but not between QCp.niab.2A and QCp.niab.4D. The two plant height QTL, QPh.niab.4B and QPh.niab.4D, which mapped to the same locations as QCp.niab.4B and QCp.niab.4D, also displayed significant genetic interactions.

A computational method for automated characterization of genetic components.

The ability to design and construct synthetic biological systems with predictable behavior could enable significant advances in medical treatment, agricultural sustainability, and bioenergy production. However, to reach a stage where such systems can be reliably designed from biological components, integrated experimental and computational techniques that enable robust component characterization are needed. In this paper we present a computational method for the automated characterization of genetic components. Our method exploits a recently developed multichannel experimental protocol and integrates bacterial growth modeling, Bayesian parameter estimation, and model selection, together with data processing steps that are amenable to automation. We implement the method within the Genetic Engineering of Cells modeling and design environment, which enables both characterization and design to be integrated within a common software framework. To demonstrate the application of the method, we quantitatively characterize a synthetic receiver device that responds to the 3-oxohexanoyl-homoserine lactone signal, across a range of experimental conditions.

Planar polarity pathway and Nance-Horan syndrome-like 1b have essential cell-autonomous functions in neuronal migration.

Components of the planar cell polarity (PCP) pathway are required for the caudal tangential migration of facial branchiomotor (FBM) neurons, but how PCP signaling regulates this migration is not understood. In a forward genetic screen, we identified a new gene, nhsl1b, required for FBM neuron migration. nhsl1b encodes a WAVE-homology domain-containing protein related to human Nance-Horan syndrome (NHS) protein and Drosophila GUK-holder (Gukh), which have been shown to interact with components of the WAVE regulatory complex that controls cytoskeletal dynamics and with the polarity protein Scribble, respectively. Nhsl1b localizes to FBM neuron membrane protrusions and interacts physically and genetically with Scrib to control FBM neuron migration. Using chimeric analysis, we show that FBM neurons have two modes of migration: one involving interactions between the neurons and their planar-polarized environment, and an alternative, collective mode involving interactions between the neurons themselves. We demonstrate that the first mode of migration requires the cell-autonomous functions of Nhsl1b and the PCP components Scrib and Vangl2 in addition to the non-autonomous functions of Scrib and Vangl2, which serve to polarize the epithelial cells in the environment of the migrating neurons. These results define a role for Nhsl1b as a neuronal effector of PCP signaling and indicate that proper FBM neuron migration is directly controlled by PCP signaling between the epithelium and the migrating neurons.

The neuroepithelial basement membrane serves as a boundary and a substrate for neuron migration in the zebrafish hindbrain.

BACKGROUND: The facial branchiomotor neurons of cranial nerve VII undergo a stereotyped tangential migration in the zebrafish hindbrain that provides an ideal system for examining the complex interactions between neurons and their environment that result in directed migration. Several studies have shown the importance of the planar cell polarity pathway in facial branchiomotor neuron migration but the role of apical-basal polarity has not been determined. Here we examine the role of the PAR-aPKC complex in forming the basal structures that guide facial branchiomotor neurons on an appropriate migratory path. RESULTS: High resolution timelapse imaging reveals that facial branchiomotor neurons begin their migration by moving slowly ventrally and posteriorly with their centrosomes oriented medially and then, upon contact with the Laminin-containing basement membrane at the rhombomere 4-rhombomere 5 boundary, speed up and reorient their centrosomes on the anterior-posterior axis. Disruption of the PAR-aPKC complex members aPKClambda, aPKCzeta, and Pard6gb results in an ectopic ventral migration in which facial branchiomotor neurons escape from the hindbrain through holes in the Laminin-containing basement membrane. Mosaic analysis reveals that the requirement for aPKC is cell-nonautonomous, indicating that it is likely required in the surrounding polarized neuroepithelium rather than in facial motor neurons themselves. Ventral facial motor neuron ectopia can be phenocopied by mutation of lamininalpha1, suggesting that it is defects in maintenance of the laminin-containing basement membrane that are the likely cause of ventral mismigration in aPKClambda+zeta double morphants. CONCLUSIONS: Our results suggest that the laminin-containing ventral basement membrane, dependent on the activity of the PAR-aPKC complex in the hindbrain neuroepithelium, is both a substrate for migration and a boundary that constrains facial branchiomotor neurons to the appropriate migratory path.

Alternative splicing of lola generates 19 transcription factors controlling axon guidance in Drosophila.

The Drosophila melanogaster transcription factor Lola (longitudinals lacking) is a pivotal regulator of neural wiring that sets the precise expression levels of proteins that execute specific axon guidance decisions. Lola has a zinc finger DNA binding domain and a BTB (for Broad-complex, Tramtrack and Bric a brac) dimerization motif. We now show that alternative splicing of the lola gene creates a family of 19 transcription factors. All lola isoforms share a common dimerization domain, but 17 have their own unique DNA-binding domains. Seven of these 17 isoforms are present in the distantly-related Dipteran Anopheles gambiae, suggesting that the properties of specific isoforms are likely to be critical to lola function. Analysis of the expression patterns of individual splice variants and of the phenotypes of mutants lacking single isoforms supports this idea and establishes that the alternative forms of lola are responsible for different functions of this gene. Thus, in this system, the alternative splicing of a key transcription factor helps to explain how a small genome encodes all the information that is necessary to specify the enormous diversity of axonal trajectories.